Lower ground state due to counter-rotating wave interaction in trapped ion system 
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We consider a single ion confined in a trap under radiation of two traveling waves of lasers. In 
the strong-excitation regime and without the restriction of Lamb-Dicke limit, the Hamiltonian of 
the system is similar to a driving Jaynes-Cummings model without rotating wave approximation 
(RWA). The approach we developed enables us to present a complete eigensolutions, which makes 
it available to compare with the solutions under the RWA. We find that, the ground state in our 
non-RWA solution is energically lower than the counterpart under the RWA. If we have the ion in the 
ground state, it is equivalent to a spin dependent force on the trapped ion. Discussion is made for 
the difference between the solutions with and without the RWA, and for the relevant experimental 
test, as well as for the possible application in quantum information processing. 

PACS numbers: 32.80.Lg, 42.50.-p, 03.67.-a 

I. INTRODUCTION 

Ultracold ions trapped as a line are considered as a promising system for quantum information processing [l[ . Since 
the first quantum gate performed in the ion trap 0, there have been a series of experiments with trapped ions to 
achieve nonclassical states [1| , simple quantum algorithm [i[ , and quantum communication Q . 

There have been also a number of proposals to employ trapped ions for quantum computing, most of which work 
only in the weak excitation regime (WER), i.e., the Rabi frequency smaller than the trap frequency. While as bigger 
Rabi frequency would lead to faster quantum gating, some proposals 0, 0, HI have aimed to achieve operations in 
the case of the Rabi frequency larger than the trap frequency, i.e., the so called strong excitation regime (SER). The 
difference of the WER from the SER is mathematically reflected in the employment of the rotating wave approximation 
(RWA), which averages out the fast oscillating terms in the interaction Hamiltonian. As the RWA is less valid with 
the larger Rabi frequency, the treatment for the SER was complicated, imcomplete and sometimes resorted to 
numerics [10|. 

In addition, the Lamb-Dicke limit strongly restricts the application of the trapped ions due to technical challenge 
and the slow quantum gating. We have noticed some ideas [ill, E2 to remove the Lamb-Dicke limit in designing 
quantum gates, which are achieved by using some complicated laser pulse sequences. 

In the present work, we investigate, from another research angle, the system mentioned above in SER and in the 
absence of the Lamb-Dicke limit. The main idea, based on an analytical approach we have developed, is to check the 
eigenvectors and the eigenenergies of such a system, with which we hope to obtain new insight into the system for 
more application. The main result in our work is a newly found ground state, energically lower than the ground state 
calculated by standard Jaynes-Cummings model. We will also present the analytical forms of the eigenvectors and 
the variance of the eigenenergies with respect to the parameters of the system, which might be used in understanding 
the time evolution of the system. 

The paper is organized as follows. In Section II we will solve the system in the absence of the RWA. Then some 
numerical results will be presented in comparison with the RWA solutions in Section III. We will discuss about the 
new results for their possible application. More extensive discussion and the conclusion are made in Section IV. Some 
analytical deduction details could be found in Appendix. 

II. THE ANALYTICAL SOLUTION OF THE SYSTEM 

As shown in Fig. 1, we consider a Raman A- type configuration, which corresponds to the actual process in NIST 
experiments. Like in [13I ]. we will employ some unitary transformations to get rid of the assumption of Lamb-Dicke 
limit and the WER. So our solution is more general than most of the previous work |14j ]. For a single trapped 
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ion experiencing two off- resonant counter-propagating traveiing wave lasers with frequencies oj\ and 0J2 , respectively, 
and in the case of a large detuning 5, we have an effective two-level system with the lasers driving the electric-dipole 
forbidden transition \g) <-> |e) by the effective laser frequency lol = uj\—u>2- So we have the dimensionless Hamiltonian 

H = ^a z + ta + ^{o+e^ + a-e-**), (1) 

in the frame rotating with u)l, where A = (loq — uj£)/v, ujq and v are the resonant frequency of the two levels of the 
ion and the trap frequency, respectively. Q is the dimensionless Rabi frequency in units of v and 77 the Lamb-Dickc 
parameter. cr±. z are usual Pauli operators, and we have x = + a for the dimensionless position operator of the ion 
with and a being operators of creation and annihilation of the phonon field, respectively. We suppose that both f2 
and v are much larger than the atomic decay rate and the phonon dissipative rate so that no dissipation is considered 
below. 

Like in [l3j . we first carry out some unitary transformations on Eq. (1) to avoid the expansion of the exponentials. 
So we have 

H 1 = UHU f = —a z + a f a + g(a) + a)a x + ea x + g 2 , (2) 



jj = J_ e W«/a ( F^rj) F{rf) \ 



where 

"~7f ' \-F^rj) F(r{))> 

with F(rf) — exp[ir](a' + a)/2], g = r//2, and e = —A/2. Eq. (2) is a typical driving Jaynes-Cummings model 
including the counter-rotating wave terms. In contrast to the usual treatments to consider the Lamb-Dicke limit 
by using the RWA in a frame rotation, we remain the counter-rotating wave interaction in the third term of the 
right-hand side of Eq. (2) in our case. To go on our treatment, we make a further rotation with V = exp (ina y /4), 
yielding 

H = VHW 1 = --cr^ + a) a + g{a) + a)o z + ea z + g 2 , (3) 

where we have used exp (i9a y )a x exp (— i9a y ) = cos(29)a x + sin(29)a z , and exp (iOa y )a z exp (— i9a y ) = cos(2#)er 2 — 
sm(29)a x . For convenience of our following treatment, we rewrite Eq. (3) to be 

H = e(| e )(e| - \g){g\) - |(|e)( 5 | + \g){e\) + Ja + gtf + a)(\e){e\ - \g)(g\)+g 2 . (4) 
Using Schrodinger equation, and the orthogonality between |e) and \g), we suppose 

\) = \<Pi)\e) + \<P2)\g), (5) 

which yields 

e\cpi) +a^a\ip 1 ) + g(a) + a)\<px) - ~|p 2 ) + g 2 \vi) = E\<fi), (6) 

- e\ip 2 ) + <Ja\ip 2 ) - g{<J + a)\tp 2 ) - —\<pi) + g 2 \tp 2 ) = E\ip 2 ). (7) 

To make the above equations concise, we apply the displacement operator D{g) = exp [g(a' — a)] on a) and a, which 
gives A = b(g)Ub{g) = a+g, = b( g y^D(g) = a^+g, B = b{-g)^aD(-g) = a-g, and fit = £)(-g)^D(-g) = 
— g. So we have 

{A^A + e)\ Vl )-^\<p 3 )=E\ Vl ), (8) 
{B^B-e)\y 2 )-^x)=E\<p 2 ). (9) 
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Obvious, the new operators work in different subspaces, which leads to different evolutions regarding different internal 
levels \g) and |e). We will later refer to this feature to be relevant to spin-dependent force. The solution of the two 
equations above can be simply set as 

N 

\<p 1 ) = J2 c n\n) A) (10) 

n=0 



I Pa) = ^2d n \n) B , (11) 

with N a large integer to be determined later, \n) A = -^{a^ + g) n \0) A = 7=f( at + g) n D(g)^\0) = ^{a^ + 
ff) w eap{- ff ot-^/2}|0), and \n) B = -^(at - g) n \0) B = ^( a t _ g)»D(-g)t\0) = ^( a t - g) n exp^at - //2}|0). 
Taking Eqs. (10) and (11) into Eqs. (8) and (9), respectively, and multiplying by a(™| and s(m|, respectively, we 
have, 

(to + e)c m - - ^2{-l) n D mn d n = Ec m , (12) 

n=0 



(to - e)d m - -^2(~l) m D mn c n = Ed mi (13) 
n=Q 

where we have set (— l) n D mn —a (m\n) B and (— l) m D mn = B (m\ri)A, whose deduction can be found in Appendix. 
Diagonizing the relevant determinants, we may have the eigenenergies Ei and the eigenvectors regarding d n and d l n 
(n = 0, • • • , N, i = 0, • • • , N). Therefore, as long as we could find a closed subspace with c l N+1 and d l N+1 approaching 
zero for a certain big integer N, we may have a complete eigensolution of the system. 



III. DISCUSSION BASED ON NUMERICS 



Before doing numerics, we first consider a treatment by involving the RWA. As the RWA solution could present 
complete eigenenergy spectra, it is interesting to make a comparison between the RWA solution and our non-RWA 
one. We consider a rotation in Eq. (2) with respect to exp{— i[{fl/2)a z + a^a]i}, which results in 

H A = + a ] a + g{aa + + a 1 a-) + g 2 , (14) 
where the RWA has been made by setting fi = 1, and we have corresponding eigenenergies 

E± = (n + g 2 + 1/2) ± gV^+T. (15) 

So the system is degenerate in the case of i] = and there are two eigenenergy spectra corresponding to as long 
as i] 7^ 0. 

Figs. 2(a) and 2(b) demonstrate two spectra, respectively, and in each figure we compare the differences between the 



RWA and non-RWA solutions |15(. In contrast to the two spectra in the RWA solution, the non-RWA solution includes 
only one spectrum. Comparing the two eigensolutions, we find that the even- number and odd- number excited levels 
in the non-RWA case correspond to E+ and E~ of the RWA case, respectively, and the difference becomes bigger and 
bigger with the increase of rj. ft is physically understandable for these differences because the RWA solution, valid only 
for small to does not work beyond the Lamb-Dicke regime. Above comparison also demonstrates the change of the ion 
trap system from an integrable case (i.e., with RWA validity) to the non-integrable case (i.e., without RWA validity). 
But besides these differences, we find an unusual result in this comparison, i.e., a new level without the counterpart 
in RWA solution appearing in our solution, which is lower than the ground state in RWA solution by v + xr\ with x 
a jy-dependent coefficient. In the viewpoint of physics, due to additional counter-rotating wave interaction involved, 
it is reasonable to have something more in our solution than the RWA case, although this does not surely lead to 
a new level lower than the previous ground state. Anyway, this is a good news for quantum information processing 
with trapped ions. As the situation in SER and beyond the Lamb-Dicke limit involves more instability, a stable 
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confinement of the ion requires a stronger trapping condition. In this sense, our solution, with the possibility to have 
the ion stay in an energically lower state, gives a hope in this respect. We will come to this point again later. 

Since no report of the new ground state had been found either theoretically or experimentally in previous pub- 
lications, we suggest to check it experimentally by resonant absorption spectrum. As shown above, in the case of 
non-zero Lamb-Dicke parameter, the degeneracy of the neighboring level spacing is released, and the bigger the 77, 
the larger the spacing difference between the neighboring levels. Therefore, an experimental test of the newly found 
ground state should be available by resonant transition between the ground and the first excited states in Fig. 2, once 
the SER is reached. We have noticed that the SER could be achieved by first cooling the ions within the Lamb-Dicke 
limit and under the WER, and then by decreasing the trap frequency by opening the trap adiabatically p. 

Since it is lower in energy than the previously recognized ground states, the new ground state we found is more 
stable, and thereby more suitable to store quantum information. Once the trapped ion is cooled down to the ground 
state in the SER, it is, as shown in Eq. (5) with n — 0, actually equivalent to the effect of a spin-dependent force on 
the trapped ion [r|. If we make Hadamard gate on the ion by \g) — ► (\g) + |e))/V2 and |e) — > (\g) — \e))/s/2, we 
reach a Schrodinger cat state, i.e., (l/2){[D+(sf)|0) + D^(-g)\0)]\g) - [D^(g)\0) - (- g)\0)]\e)} . Two ions confined 
in a trap in above situation will yield two-qubit gates without really exciting the vibrational mode . It is also the 
way with this spin-dependent force towards scalable quantum information processing fl2| . As in SER, we may have 
larger Rabi frequency than in WER, the quantum gate could be in principle carried out faster in the SER. 

In addition, as it is convergent throughout the parameter subspace, our complete eigensolution enables us to 
accurately write down the state of the system at an arbitrary evolution time, provided that we have known the initial 
state. This would be useful for future experiments in preparing non-classical states and in designing any desired 
quantum gates with trapped ions in the SER and beyond the Lamb-Dicke limit. Moreover, as shown in Figs 3(a), 
3(b ) and 3(c), our present solution is helpful for us to understand the particular solutions in previous publication 
The comparison in the figures shows that the results in [l3| are actually mixtures of different eigensolutions. For 
example, the lowest level in Fig. 2 in [TH, corresponding to f2 = 2 and ry = 0.2, is actually constituted at least by the 
third, the fourth, and the fifth excited states of the eigensolution. 

IV. FURTHER DISCUSSION AND CONCLUSION 

The observation of the counter-rotating effects is an interesting topic discussed previously. In [ijj , a standard method 
is used to study the observable effects regardin g th e rotating and the counter-rotating terms in the Jaynes-Cummings 
model, including to observe Bloch-Siegert shift [181 ] and quantum chaos in a cavity QED by using differently polarized 
lights. A recent work [19( for a two-photon Jaynes-Cummings model has also investigated the observability of the 
counter-rotating terms. By using perturbation theory, the authors claimed that the counter-rotating effects, although 
very small, can be in principle observed by measuring the energy of the atom going through the cavity. Actually, for 
the cavity QED system without any external source involved, it is generally thought that the counter-rotating terms 
only make contribution in some virtual fluctuations of the energy in the weak coupling regime. While the interference 
between the rotating and counter-rotating contributions could result in some phase dependent effects 20]. Anyway, 
if there is an external source, for example, the laser radiating a trapped ultracold ion, the counter-rotating terms will 
show their effects, e.g., related to heating in the case of WER [2l|. In this sense, our result is somewhat amazing 
because the counter-rotating interaction in the SER, making entanglement between internal and vibrational states of 
the trapped ion, plays positive role in the ion trapping. 

We argue that our approach is applicable to different physical processes involving counter-rotating interaction. Since 
the counter-rotating terms result in energy nonconservation in single quanta processes, usual techniques cannot solve 
the Hamiltonian with eigenstates spanning in an open form. In this case, path-integral approach [22[ and perturbation 
approach (2p| , assisted by numerical techniques were employed in the weak coupling regime of the Jaynes-Cummings 
model. In contrast, our method, based on the diagonalization of the coherent-state subspace, could in principle study 
the Jaynes-Cummings model without the RWA in any cases. We have also noticed a recent publication [23| to treat 
a strongly coupled two-level system to a quntum oscillator under an adiabatic approximation, in which something is 
similar to our work in the solution of the Hamiltonian in the absence of the RWA. But due to the different features 
in their system from our atomic case, the two-level splitting term, much smaller compared to other terms, can be 
taken as a perturbation. So the advantage of that treatment is the possibility to analytically obtain good approximate 
solutions. In contrast, not any approximation is used in our solution, which should be more efficient to do the relevant 
job. 

In summary, we have investigated the eigensolution of the system with a single trapped ion, experiencing two 
traveling waves of lasers, in the SER and in the absence of the Lamb-Dicke limit. We have found the ground state 
in the non-RWA case to be energically lower than the counterpart of the solution with RWA, which would be useful 
for quantum information storage and for quantum computing. The analytical forms of the eigenfunction and the 
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complete set of the eigensolutions would be helpful for us to understand a trapped ion in the SER and with a large 
Lamb-Dicke parameter. We argue that our work would be applied to different systems in dealing with strong coupling 
problems. 
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VI. APPENDIX 

We give the deduction of A{m\n)B and B{m\n)A below, 

A (m\n) B = — L=(0| e -^-s 2 /2 (a + ff )™( a t _ ff )«X-s 2 /2| 0) 
\Jrti\n\ 

= -jl=e-^ 2 (0\(a + g) m e^e-^(a^ - g) n \0) 
Vmini 

= - 7 I= e - 2s2 (0|(a + 2. 9 r(at -2g) n \0) = (-l) n D mn , 

with 

min[m,n] , j— - \ m+n -2i 

mn " ^ 1 j {m-i)\{n-i)\i\ ' 

It is easily proven following a similar step to above that 

B (m\n) A = ^^( | e ^ 2 /2 (a _ gr(a t + gre - 9 at- 9 2 /2| 0); 
Vm'.ni 

would finally get to (— l) m D mn . 
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The captions of the figures 

Fig. 1 Schematic of a single trapped ion under radiation of two traveling wave lasers, where lo\ and 0J2 are frequencies 
regarding the two lasers, respectively, loq is the resonant frequency between \g) and |e), and S and A are relevant 
detunings. This is a typical Raman process employed in NIST experiments, with for example Be + , for quantum 
computing. 

Fig. 2 The eigenenergy spectra with f2 = 1, where (a) and (b) correspond to two different sets of eigenenergies 
with respect to Lamb-Dicke parameter. In (a) the comparison is made between £"+ in the RWA case (dashed-dotted 
curves) and E n with n = even numbers in the non-RWA case (star curves for n = and solid curves for others); In 
(b) the comparison is for E~ in the RWA case (dashed-dotted curves) to E n with n = odd numbers in the non-RWA 
case (solid curves). 

Fig. 3 The eigenenergy with respect to the detuning A, where for convenience of comparison we have used the same 
parameter numbers as in (l3j . For clarity, we plot the different levels with different lines. The parameter numbers are 
ft = 2, and (a) 77 = 0.2; (b) 77 = 0.4; (c) 77 = 0.6. 



